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A great deal of work has been published on the capillarity-driven motion of an 
ordinary liquid over a solid surface (the isothermal dynamic contact line problem). The 
solution to this problem may be stated in terms of a functional relationship between the 
contact line speed U, the dynamic contact angle 0, and, the relevant properties of the 
system's three phases and their interfaces. One such relationship is the much-quoted 
empirical correlation of Hoffman (1975), who showed that for Ca=|iU/a<l, where |i is 
the liquid viscosity and o the surface tension, a large amount of data spanning many 
orders of magnitude in Ca could be correlated approximately by the formula Ca“k(0 3 - 
0 e 3 ), where 0 e is the equilibrium contact angle and k=0.013. On the theoretical side, the 
continuum theory of the moving contact line contains the well-known non-integrable 
shear stress singularity at the contact line [e.g. Tanner (1979), de Gennes (1985, 1990), 

Cox (1986), Leger (1992)]. The singularity can be removed by specifying a velocity slip 
model or cut-off very close to the contact line [e.g. Huh (1971, 1977), Dussan V. (1974, 
1976, 1979, 1991), Durbin (1988), Haley (1991), Hocking (1992)], and the parameters 
can be chosen so as to yield agreement with experiment [Fermigier (1991), Chen (1992), 
Leger (1992)]. An entirely different theory for the dynamic contact line has been put 
forward by Blake (1969) based on Eyring's reaction rate theory for liquids. With an 
appropriate selection of physical parameters, this molecular-kinetic theory can also be 
brought into agreement with experimental data [Blake (1969, 1979, 1994), Gribanova 
(1992)]. Some progress has also been made recently in applying molecular-dynamic 
computations to the contact line (Thompson et al, 1993; Koplic & Banavar, 1995), but 
these are limited to mesoscale systems. Clearly, despite the large amount of work, this 
problem is not yet conclusively resolved. 

In contrast to the ordinary contact line problem, virtually no information is 
available on the similar problem associated with a molten material spreading on a solid 
which is below the melt's fusion point. The latter is a more complex problem in which 
heat transfer and solidification take place simultaneously with spreading, and requires 
answers not only for the hot melt’s advance speed over the cold solid as a function of 
contact angle, but also for how one is to predict the point of the molten contact line's 
arrest by freezing. These issues are of importance in evolving methods of materials 
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processing. For example, finely featured components can be fabricated under computer 
control by precise deposition of individual molten microdrops, each droplet solidifying 
by heat loss after it lands (Gao & Sonin, 1994). In this case, the point of a droplet's 
contact line arrest will control its post-deposition solid shape, that is, the shape of the 
basic microscale "building block" of the process. 

The purpose of our work is to develop, based on both experiments and theory, an 
understanding of the dynamic processes that occur when a molten droplet touches a 
subcooled solid, spreads partly over it by capillary action, and freezes. We seek answers 
to the following basic questions. First, what is the relationship between the melt’s contact 
line speed and the apparent (dynamic) contact angle? Secondly, at what point will the 
contact line motion be arrested by freezing? 

A simple well-defined example is the deposition of a molten droplet at 
temperature T 0 onto a semi-infinite, solid target at a temperature T t below the melt's 
freezing point Tf (Fig. 1). If the scale is so small, or gravity so weak, that gravitational 
effects are negligible (small Bond number), and the deposition so gentle that the molten 
droplet does not break up upon impact (small or moderate Weber number), the droplet 
will typically spread in liquid form until the contact line is arrested by freezing, after 
which point the droplet's footprint will remain fixed, though the liquid above it may still 
undergo dynamic oscillations until solidification is complete. Since the subcooled target 
has an infinite thermal capacity, the final state will always be one where the deposited 
material is solid and at the target temperature. 

The talk will describe three components of our work: (i) deposition experiments 
with small molten droplets, (ii) investigation of the dynamics of the molten contact line 
by means of a novel forced spreading method, and (iii) an attempt to provide a theoretical 
framework for answering the basic questions posed above. 

(i) Deposition experiments were performed with molten droplets of waxes, water, 
and mercury. The thermal-dynamic scaling parameters for deposition are listed in Fig. 2, 
simplified to the case where the Weber number is small and the molten droplet and solid 
target are of the same material. Fig. 3 gives material properties and parameter values. 
The droplet radii ranged from 25 pm for the wax to about 300 pm for mercury and up to 
1000 pm for water. Mercury was included because it possesses typical metal properties 
but, owing to its low melting point, is relatively free of interfacial oxidation problems. A 
high speed camera provided information on the global droplet dynamics during impact 
and spreading, including time scales and the final drop shape (Fig. 4 gives the ratio of the 
arrested footprint radius R* and a, from which contact angle can be derived), but did not 
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give enough resolution to accurately determine dynamic contact angles during the 
spreading process. 

(ii) To better investigate the contact line dynamics, we devised the forced 
spreading experiment shown schematically in Fig. 5. Using an apparatus similar to Gao 
& Sonin’s (1994), a continuous stream of molten microdrops (radius 25 p.m) is ejected at 
frequency f and ballistically delivered to a point on the target surface. The droplets 
coalesce as they arrive to form a sessile liquid macrodrop which grows and spreads 
laterally over the target. This simulates the spreading of a single, impacting droplet. The 
delivery frequency provides control over the contact line velocity U=dR/dt, where R is 
the macrodrop's contact line radius. The contact line motion can be filmed with a 
standard CCD camera through a microscope objective. 

The video sequence of Fig. 6 shows a forced spreading experiment with molten 
wax on solid wax (only waxes have been studied with this technique so far). In the early 
stages the contact line moves much as it would for an ordinary liquid, and the relationship 
between 0 and U is well represented by Hoffman's law with 0 C =O. The contact line speed 
decreases as the macrodrop's base radius grows and, after the advance velocity slows to a 
critical value U*, the contact line arrests rather abruptly at a radius R* (Fig. 7). Our 
experiments show that U* depends on target subcooling (Fig. 8), but is essentially 
independent of melt superheat (T 0 -Tf)/(Tf-Tt) and independent also of the macrodrop's 
spreading history (Fig. 8), which in this experiment is controlled via the delivery rate f 
and time. The primary parameter which determines U* appears to be the Stefan number 
S=c(Tf-Tt)/L based on target subcooling, c being the specific heat and L the latent heat. 
Contact line arrest is thus controlled mainly by processes near the contact line. 

(iii) A freezing mechanism may be postulated based on the angle 0 S of the 
solidification front at the contact line (Fig. 9). A melt droplet will spread until the 
apparent liquid contact angle 0 approaches 0 S , at which point the supply of melt to the 
contact line cannot be maintained and the spreading must stop. If this is correct, the 
melt's apparent contact angle 0 will be given approximately by Hoffman's law with 0 e =O 
as long as 0s«0, and contact line freezing will occur when 0 S approaches 0. The arrest 
condition would therefore be obtained by calculating 0 S . 

There is one problem with this: if we use the conventional continuum equations 
and boundary conditions (Fig. 10), the heat flux which controls the angle of the fusion 
interface is singular at the contact line (it has the same power-law form near the CL as 
Anderson & Davis's (1994) analysis for the purely conductive case). The angle 0 S cannot 
therefore be computed unless the continuum equations and/or boundary conditions are 
modified with a physically correct mechanism that removes the mathematical singularity. 
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The mechanism is at this time open to some speculation. Our first approach has been to 
ask whether a material can be characterized with a given cutoff length X. We ask, at what 
distance from the contact line should the heat flux of the continuum solution be 
terminated to give a 0 S value in agreement with the experimental data for the arrest 
condition? Computations (Figs 10 and 11) show that the data for each of the three 
materi als we have studied can indeed be adequately matched with a cut-off length X, the 
value of which ranges from 0.08 to 0.3 pm depending on material (Fig. 12). These 
lengths are large compared with molecular scale, but smaller than the typical velocity slip 
lengths that have been used to fit continuum theories of the ordinary dynamic contact line 
to experimental data. 

This research is at an early stage, and further work is required on both the 
experimental and theoretical sides to resolve the heat flux singularity and lay the basis for 
a complete theory of the molten contact line. 

Acknowledgement: This work was supported by the National Science Foundation under 
Grants CTS-9122123 and CTS-9413026. 
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Fig. 1: Molten droplet deposition on subcooled solid of its 
own kind 
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Fig.2: Global scaling parameters for capillarity- 
driven droplet deposition 
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Fig. 3: Properties of Selected Materials 
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Fig. 4: Effects of thermal conditions on final droplet shape 



• decreases with increasing target subcooling 
S ■ c(Tf-Ti) / L 

• is insensitive to melt superheat p * (T 0 -Tf) / (Tf-Tt) 
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Fig. 5: The forced spreading experiment 
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Fig. 6: Video sequence of forced spreading 
experiment. Microcrystalline wax, T 0 =117 r C, 
T t =80 X, Tf=91 X, f= 1kHz, the rectangle's width 
is 100 [im 



Fig. 7: Results from forced spreading 
experiments 
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Contact line arrest conditions (apparent contact angle 6s 
and speed U* just prior to arrest) are: 

(i) independent of macrodrop spreading history (i.e. 
no dependence on frequency of deposition) 

(ii) dependent on target subcooling 



Fig. 8: CL speed just prior to arrest increases 
with target subcooling S and is independent of 
feed frequency f 
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Fig. 9: The solidification front at the 
contact line 


Fig. 10: Calculation of heat flux near the CL 
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Fig. 11: Finite element calculation of temperature Fig. 12: Comparison of experiments 

field for a wedge domain with 0 = n/6 and theory on molten contact line arrest 







